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R. Maiolino^, a. Alonso-Herrero^, S. Anders*, A. Quillen^, M.J. Rieke^, G.H. Rieke^, 

L.E. Tacconi-Garman" 






O 

00 



ABSTRACT 

We report the discovery of gas inflow motions towards the active nucleus of the Circinus galaxy caused 
by the non-axisymmetric potential of a nuclear gas bar. Evidence for dust associated with the bar comes 
from the HST/NICMOS H — K color map, whereas the streaming motions along the gas bar are seen 
in the velocity field of the H2S(l)(l-0) emission line. The gas bar is about lOOpc long with a visual 
extinction in excess of 10 mag. Indication for the gaseous nature of this bar comes from the lack of a 
stellar counterpart even in the K band where the extinction is greatly reduced. 

We also use the NICMOS emission line images (Paa, [Sivi], and [Fell]) to study the innermost region 
of the ionization cones and the nuclear star forming activity. We discuss the possible relationship of 
these components with the gaseous bar. 

Subject headings: Galaxies: individual: Circinus ~ galaxies: nuclei ~ galaxies: active - galaxies: 
Seyfert - galaxies: kinematics and dynamics - infrared: galaxies 
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1. INTRODUCTION 

Active galactic nuclei (AGNs) are thought to be powered 
by an accreting supermassive black hole (Mbh ~ 10® — 
10^°Mq). The required accretion rate ranges from lO^'^ — 
10-1 Moyr-i (Seyfert nuclei) up to 10 - 100 Mgyr^i 
(QSOs). Although the gas content of host galaxies is 
generally enough to fuel AGNs over the Hubble time, the 
mechanism responsible for transporting gas from the kpc 
scale into the nuclear region is still unclear. Indeed, the 
galactic gas must lose about 99.9% of its angular momen- 
tum before reaching the central few parsecs, where vis- 
cosity can be effective. The required accretion rates are 
lower in Seyfert galaxies than in QSOs and the gas inflow 
requirements from the host galaxy are relaxed. However, 
because of the large number of Seyfert galaxies (~ 10%) 
in the local Universe and the even larger fraction of AGNs 
at high redshifts (Ho et al. 1997, Maiohno & Rieke 1995, 
Hasinger et al. 1998, Miyaji et al. 1999), some mechanism 
for transporting gas from the circumnuclear 100 pc into 
the central parsec is still required to maintain the activity 
of these nuclei. Also, the accumulation of large quantities 
of gas in the innermost region of Seyfert nuclei is inferred 
from the observation that most of them are obscured by 
gaseous column densities in excess of 10^''cm~^ within the 
central 10 pc (Risaliti et al. 1999) and from the detection 
of maser emission within the nuclear few parsecs (Green- 
hiU et al. 1997a, 1997b, 1998b, GreenhiU & Gwinn 1997, 
Miyoshi et al. 1995). 

Non-axisymmetric gravitational potentials are thought 
to be an efficient mechanism for transporting gas from the 
host galaxy into the central region (Athanassoula 1992, 
Heifer & Blitz 1995, Laine et al. 1999, Hernquist 1989, 
Quillen et al. 1995, Regan et al. 1997, Shlosman et al. 
1990, Maiohno et al. 1999a). However, within the cen- 



tral few hundred parsecs large scale disturbances have lit- 
tle effect in removing angular momentum from the gas, 
since on these scales the gravitational field is dominated 
by the central region of the bulge that is axially symmet- 
ric. Nested, secondary stellar bars have been observed in 
several barred system (Jungwiert et al. 1997, Friedli ct 
al. 1996, Wozniak et al. 1995). These secondary bars are 
relatively small (a few 100 pc long) and could be respon- 
sible for transporting gas into the nuclear region to feed 
a black hole. Indeed, the fraction of Seyfert nuclei among 
the 40 double-barred systems known so far is about 1/3 to 
be compared to the fraction of Seyferts among all galaxies 
that is about 1/10 (Friedli 1999). Alternatively, Shlosman 
et al. (1989), Wada & Habe (1992) and Heller & Shlosman 
(1994) proposed that if the mass of gas gathered into the 
central region exceeds about 10% of the dynamical mass 
then the gaseous disk develops gravitational instabilities 
and forms a gaseous bar, which can drive gas further into 
the nuclear region. Indeed, there is evidence in several 
Seyfert nuclei for concentrations of gas whose mass ap- 
proaches the dynamical mass (Meixner et al. 1990, Henkel 
et al. 1991, Risaliti et al. 1999), therefore gaseous bars 
are expected to form and to provide a mechanisms for the 
fuelling of the nuclear region. However, neither of these 
mechanisms appear to operate in the majority of Seyfert 
galaxies. 

Understanding the mechanism responsible for the fu- 
elling of active galaxies is also relevant to other issues. The 
dense gas tunneled into the central region is photoionized 
by UV radiation from the nuclear source and/or is shock 
excited by the interaction with radio jets. Also, the large 
mass of gas driven into the central region and the perturba- 
tions due to the fuelling mechanism are expected to boost 
the star forming activity in the central region. There- 
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fore, the fuelling mechanism of AGNs is likely to play a 
key role in the connection between starbursts and AGNs 
(Gonzalez-Dclgado et al. 1998, Maiolino et al. 1997). 

The Circinus galaxy is a nearby (4 Mpc), edge-on 
(incl.^ 65°), Sb~d system that is seen through a low in- 
terstellar extinction window near the Galactic plane (Ay 
= 1.5 mag, Freeman et al. 1977). The nuclear optical 
line ratios (Oliva et al. 1994) are typical of a Seyfert 2 
galaxy. This classification is supported by the detection of 
intense coronal lines (Oliva et al. 1994, Moorwood et al. 
1996), the discovery of an intense X-ray iron 6.4 keV line 
(Matt et al. 1996), rapid variation of the powerful H2O 
maser emission (Greenhill et al. 1997c) and a prominent 
ionization cone in the [OIIIJA5007 maps (Marconi et al. 
1994a) with filamentary supersonic outflows (Veilleux & 
Bland-Hawthorn 1997). The detection of an X-ray excess 
above 30 keV reveals that the central engine is heavily ob- 
scured by a gaseous column density Nh = 4 x lO^'^cm^^ 
along our line of sight (Matt et al. 1999). Large amounts 
of molecular gas in the central region of this galaxy have 
been inferred also from CO observations (Aalto et al. 1995, 
Elmouttie et al. 1998a, Curran et al. 1998, Johansson et 
al. 1991) and from the detection of a nuclear maser disk 
(Greenhill et al. 1998). Ha and [SII] narrow band images 
(Marconi et al. 1994a) have revealed a star forming ring at 
a radius of -10"=200pc, while Maiohno et al. (1998) have 
found evidence for a young stellar population in the cen- 
tral 10-100 pc with typical ages of the order of 10^ yr. As 
a consequence, Circinus is an excellent candidate to search 
for a gaseous bar that might be feeding both the black hole 
and the star forming activity in the central region. 

In this paper we present high angular resolution near in- 
frared HST images and integral field spectra of the central 
region of the Circinus galaxy that reveal the presence of a 
small nuclear gaseous bar and show evidence for streaming 
motions of the gas along the bar. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. NICMOS observations 

ff5'T/NICM0S observations of Circinus were obtained 
during two different orbits on March 16 and October 16 
1998 using the NIC2 and NIC3 cameras. In Table 1 we 
give the log of the NICMOS observations. Column (1) lists 
the NICMOS camera; column (2), the filter; column (3), 
the corresponding ground-based broad-band filter, emis- 
sion line or adjacent continuum; column (4), the total in- 
tegration time in seconds; and column (5), the original 
orientation of the images. The observations were taken 
in a spiral dither with a 5.5 pixel spacing between each 
of four positions. The plate scales for cameras NIC2 and 
NIC3 are 0.076" /pixel and 0.204" /pixel respectively 

The reduction of the NICMOS images used routines 
from the package NicRed (McLeod 1997). A master dark 
image was produced by combining between 10 and 20 
darks for a given sample sequence after the subtraction 
of the first readout. The darks were taken from other 
programs carried out close in time. When possible we 
also generated our own flatfield images from on-orbit data. 
The flatfields used for the NIC3 narrow-band images were 
obtained as part of the camera 3 campaign in January 
1998 and were kindly reduced and provided by Dr. Rodger 
Thompson. Only for the NIC2 F187N and NIC2 F190N 



filters were in-flight flats not available, so thermal-vacuum 
flatflelds were used instead. The data reduction was per- 
formed using the following steps: subtraction of the first 
readout; dark current subtraction on a read-by-read basis; 
correction for linearity and cosmic ray rejection (using the 
fullfit routine in NicRed); and flatflelding. The individ- 
ual dithered galaxy images were registered to a common 
position using fractional pixel offsets and a cubic spline 
interpolation, and combined to produce the flnal image of 
each field. Since our NICMOS images were obtained af- 
ter August 1997, no correction for the pedestal effect was 
necessary. 

Prior to flux calibrating the images, the background 
needs to be subtracted from the images (this is important 
only for filters at wavelengths longer than 2^m). Due to 
the large projected size of Circinus, the field of view of the 
NIC2 images is not large enough to allow measurements 
of the background on blank corners of the images, so we 
used background measurements taken during the Servic- 
ing Mission Observatory Verification (SMOV) program for 
the filter NIC2 F222M, and from our observations of the 
interacting galaxy Arp 299 for the NIC2 F160W filter. The 
NIC3 mosaics cover a larger field of view, and an estimate 
of the thermal background was obtained from the corners 
of the images. The flux calibration of the broad-band, on- 
line, and off-line images was performed using conversion 
factors (from ADU/s to mJy) from measurements of the 
standard star P330-E during SMOV. 

NICMOS only provides continuum bands to the red of 
the emission lines. The extinction in the nuclear regions 
of Circinus is very high (even at infrared wavelengths), 
and a straight subtraction of the longer wavelength off- 
line image results in an overcorrection of the continuum 
at the wavelength of the emission line. To estimate the 
continuum at 1.644 ^m, 1.96 /im, 2.122 /im and 1.87 /xm 
for the NIC3 F164N, NIC3 F196N and NIC3 212N respec- 
tively, we fit the continuum between 1.66 /im and 1.90 /xm 
using the NIC3 F166N and NIC3 F200N fine-free images 
for the first two, and the NIC3 F200N and NIC3 F215N 
for the last image. The same procedure was followed to 
subtract the contimuim at 1.87 /im for the NIC2 F187N 
image, using the NIC2 F160W and F222M images. The 
continuum fit is simply a linear regression of the flux as 
a function of the wavelength for each point in the im- 
age. An image at the required continuum wavelength is 
then constructed by extrapolation or interpolation of the 
corresponding fit. The extrapolated continuum is sub- 
tracted from the line+continuum images to produce the fi- 
nal [Fen] 1.644 /im, Paa, [Si vi]1.96/im and H2 at 2.121 /im 
images. Due to the low surface brightness of the H2 emis- 
sion together with problems associated with the continuum 
subtraction, the resulting H2 image had relatively low sig- 
nal to noise, and therefore we used the H2 map obtained 
with the 3D integral field spectrometer (see next section). 
Note however that when we rebinned the NICMOS H2 im- 
age to the angular resolution of the 3D image, both images 
showed a similar morphology. 

The images presented in this paper were taken during 
different orbits and hence have different orientations (see 
Table 1). Image rotation was removed by linear interpola- 
tion to the conventional orientation (north up, east to the 
left), both to provide the same orientation for all the im- 
ages, and also to allow easy comparisons with previously 
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Table 1 
Log of the NICMOS observations. 



Camera 



Filter 



line/broad/cont 



in 



orient 



NIC2 
NIC2 
NIC2 
NIC2 

NIC3 
NIC3 
NIC3 
NIC3 
NIC3 
NIC3 



F160W 
F222M 
F187N 
F190N 
F164N 
F166N 
F196N 
F200N 
F212N 
F215N 



H 

K 

Paa 

Paa cont. 

[Fell] 1.664 /im 

[Fell] cont. 

[Si VI] 1.96 /im 

[SiVl] cont. 

H2 2.122 ^m 

H2 cont. 



192 
192 
384 
384 
256 
256 
320 
320 
640 
640 



94° 

94° 

-64° 

-64° 

94° 

94° 

94° 

94° 

-64° 

-64° 



published data. 

2.2. 3D integral field spectroscopy 

We observed the nuclear region of the Circinus galaxy 
on March 18"^ and 19'^ 1998 with the near infrared in- 
tegral field spectrometer 3D (Weitzel et al. 1996), as- 
sisted by ROGUE, a first order adaptive optics system 
(Thatte et al. 1995), at the Anglo Australian Telescope. 
3D slices the focal plane into 16 slits and disperses their 
light in wavelength; the spectrum of the slices is then im- 
aged onto a NICM0S3 detector (256x256 pixels). The 
resulting "cube" provides simultaneous spectra of an area 
6.4" X 6.4", as projected on the sky, divided into 16x 16 pix- 
els (i.e. 0.4"/pix). Circinus was observed with a grism pro- 
viding a spectral resolution of A/AA = 2100 and covering 
the spectral range from 1.95/j,m to 2.21/j,m. The average 
seeing for accepted observations, after tip-tilt correction, 
was about 0".7 — 0".8. The (limited) field of view was not 
centered on the nucleus but was moved about 1 arcsec to 
the south-east so that it could include a dusty feature we 
were interested in. The total on source integration time 
was 20 minutes, split in twelve elementary integrations of 
100 seconds. The elementary on-source integrations were 
interleaved with an equal number of acquisitions of sky 
frames for background subtraction. To avoid and control 
instrumental artifact we dithered the field of view around 
the region of interest so that each region of the source 
was observed by different regions of the slicer and of the 
detector. We observed a nearby 09 star to correct the at- 
mospheric transmission features. Data reduction was per- 
formed as described in Weitzel et al. (1996) and Thatte et 
al. (1997). The subtraction of the continuum underneath 
the emission lines was performed with a first order fit to 
the continuum approximately ± 0.03 /im from each line. 

3. THE NUCLEAR GAS BAR 

3.1. Evidence for a radial dust lane 

Fig. la shows a false color image of the central region of 
the Circinus galaxy (26" x 29" = 530 x 590 pc) obtained 
by combining HST images at 0.606 pun (WFPC2-F606W, 
blue), at 1.66 ^m (NIC3-F166N, green), and at 2.0 /zm 
(NIC3-F200N, red). The galactic disk is inclined by about 
65° with respect to our line of sight (the South-East is the 



near side of the disk) and the position angle of the major 
axis is about 25°. The image shows a spiral pattern with 
prominent dust lanes, consistent with the large amount 
of gas in this region as inferred from the CO millimet- 
ric data. Fig. lb shows the H-K color map constructed 
with the NIC2-F160W (- H band) and NIC2-F222M (~ 
K band) images. 

The H-K color of the nucleus (identified with the 2/zm 
peak) is extremely red and in the F222M (^ K band) 
image the nucleus dominates the emission within the cen- 
tral two arcseconds. As discussed in detail in Maiolino et 
al. (1998), the powerful K band emission is due to hot 
dust, heated by the active nucleus and located in a re- 
gion smaller than 1.5 pc {— 0".07 arcsec) in radius. In 
Fig. lb we have subtracted a nuclear unresolved source 
using the NIC2 point spread function as modeled by the 
TinyTim software (Krist & Hook 1997). This subtrac- 
tion allows study of the H-K color of the circumnuclear 
two arcseconds (although in the innermost arcsecond there 
are some residuals probably due to imperfect subtraction). 
The PSF subtraction does not affect the color map outside 
this inner region. 

The H-K color of the circumnuclear region is still sig- 
nificantly redder than colors typical of stellar populations 
older than 10^ yr: (H-K)« 0.22, as observed both in galax- 
ies that do not suffer heavy extinction, and also as pre- 
dicted by population models. Such off-nuclear red colors 
cannot be ascribed to hot dust emission in the K band, 
for the CO stellar bands at 2.29 /xm do not appear diluted 
by non-stellar emission in these regions (Maiolino et al. 
1998). In addition, the H-K color does not correspond 
with the Paa and Br7 emission (see Sect.||), as expected 
if the UV continuum were heating dust locally. The only 
plausible explanation is that the H-K color excess is due 
to reddening. 

There is a moderate large-scale extinction gradient 
across the minor axis of the galactic disk (i.e. from the 
NW to the SE), resulting from the greater effectiveness of 
dust absorption on the near side of the galactic disk (fewer 
bulge stars in front of us and more bulge stars behind the 
dusty disk). However, the most prominent feature of the 
"extinction map" in Fig. lb is the L-shaped feature ex- 
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Fig. 1. — a) True color image of the central region of the Circinus galaxy obtained by combining HST images at 0.606 fim (WFPC2-F606W, 
blue), at 1.66 fim (NIC3-F166N, green) and at 2.0 ^m (NIC3-F200N, red), b) H-K color map of the nuclear region constructed with the 
NIC2-F160W (~ H band) and NIC2-F222M (^ K band) images. The nuclear PSF was subtracted in both filters before producing the color 
map. The right hand side bar converts the color code into H-K color excess and into visual extinction in the case of a foreground screen 
model. 
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Fig. 2. — Schematic drawing of the geometry of the various components within the central few 100 pc of Circinus. 
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Fig. 3. — R-H color map obtained by means of the WFPC2-F606W image and the NIC2-F160W image. The straight lines indicate the 
orientation of the light cone observed in [OIII]. 



tending from the nucleus to the south-east^: we beheve 
that the section of this feature elongated radially to the 
south is tracing the near side of a nuclear gaseous-dusty 
bar. The linear projected (half) size is about 50 pc. The 
southern end of the dusty bar tilts by about 90 degrees 
to the east, in a tightly wound spiral arm-like feature (see 
also Fig. la); this morphology is similar to that observed 
in stellar bars on larger scales, but in a gaseous version. 
A stellar counterpart is not present even in the K band 
image, where the extinction is reduced, suggesting the bar 
is purely gaseous. 

An alternative explanation of the radial dust lane could 
be a warped disk (eg. Quillen et al. 1993). However, 
as discussed in the next section, the velocity field of the 
molecular gas strongly deviates from circular motion, rul- 
ing out this interpretation. 

The northern extremity of the gaseous bar ends in a 
dusty feature that resembles a circumnuclear disk, about 
20 pc in radius. Only the near side of this disk-like fea- 
ture is clearly delineated, possibly as a consequence its 
being seen in absorption (as for the bar) and therefore the 
contrast is expected to be much higher on the near side. 
According to models, the gaseous bar is expected to form 
a gaseous disk in its central part during the latest stages of 
its evolution as a consequence of cloud-cloud viscous col- 
lisions at the denser regions near the aphelia of the orbits 
(Wada & Habe 1992). 

Fig. 2 schematically shows the geometry of the gaseous 
features discussed in this section along with other com- 
ponents discussed in the following sections. The ob- 
served angle between the line of the minor axis and the 
gaseous bar is aobs ~ 30°, while its deprojected value is 
a = tan~^(tanaobs cos i) — 16°, where i is the inclination 
angle of the galaxy (65°). The deprojected length of the 



bar is 



L(bar) = Lobs (bar)- 



cos aobs 



105 pc 



(1) 



cos I ■ cos a 

We can use the H-K color excess to derive the extinc- 
tion affecting the nuclear region and, in particular, the 
extinction introduced by the dust in the gas bar. By as- 
suming that the intrinsic color of the stellar population 
is (H-K)o =0.22, we derive the magnitude of the dust 
absorption affecting the stellar light. By adopting the ex- 
tinction curve of He et al. (1995) and a foreground dusty 
screen model we obtain: Eh-k = 0.065 Ay. The color bar 
on the right-hand side of Fig. lb converts the color coding 
of the H-K map into color excess Eh-k and into visual 
extinction according to the extreme case of a foreground 
dusty screen. 

The foreground screen model is the most conservative 
assumption: given a color excess it provides the lowest 
Ay. In a more realistic case, dust will be partly mixed 
with stars. Also, a fraction of stars will also be present in 
foreground and their contribution is more important in the 
most obscured regions. The latter effect is dominant in the 
optical, indeed the heavy dust obscuration affecting the 
nuclear region is so high that these regions are essentially 
invisible at optical wavelengths. This effect becomes more 
obvious in Fig. 3 where we have constructed an R-H color 
map of the central region of Circinus with the WFPC2 
F606W and the NIC2 F160W images. This map shows 
a very irregular structure and no evidence for the dusty 
bar detected in the H-K color map. This is a consequence 
of the heavy obscuration that completely absorbs the op- 
tical light and that, therefore, makes the R band image 
completely dominated by the foreground stars and more 
affected by small variations in the extinction. We tried 
to describe these effects with a simple but more quantita- 



^Evidence for this strucure, though at lower angular resolution, is observed also in the H— K color map presented in Maiolino et al. (1998). 
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Fig. 4. — Expected H— K and R-H color excess in the case of a dusty screen that obscures only a fraction of the stars (the whole system is 
then obscured by a Galactic extinction Ay = 1.5 mag). Each curve gives the color excess for an increasing absorption of the screen and for 
a given fraction of stars behind the dusty screen. Points along each curve are separated by A Ay = 1, and the total extinction ranges from 
Ay = to Ay = 20. The two boxes show the range of color excess along the gas bar and outside the bar. 



tive model. We make the simplistic assumption that the 
distribution of dust and stars can be described by a homo- 
geneous dusty screen that hides most of the stars while a 
small fraction of the stars is in the foreground. This sim- 
ple model was used for other heavily obscured systems (eg. 
McLeod et al. 1993). The whole system is then obscured 
by a Galactic absorption Ay (Gal) = 1.5 mag (Freeman 
et al. 1977). In this model, the color excess between two 
wavebands at Al and at A 2 is given by 



Eai- 



A2 



AAi(Gal) 
- 2.5^05 



AA2(Gal) 

-fhid(l-10'^ 



l-fhid(l-10-*A2(='=r)/2-5) 



(2) 



where fhid is the fraction of stars that arc hidden by the 
dust screen, while A(Gal) and A(scr) are the Galactic and 
screen absorption respectively. Fig. 4 shows the expected 
H-K and R-H color excess as a function of Av(scr) and 
for different fraction of stars behind the dusty screen. Two 
boxes indicate the color excesses measured along and out- 
side the gas bar. The R-H color excess is the more uncer- 
tain quantity, not only bacause of the observed large dis- 
persion (Fig. 3), but also because the unabsorbed, intrinsic 
R-H color is not well defined. We adopted (R — H)o — 2.3 
from an average of the colors of Sc spirals in Frogel (1985) 
and in Fukugita et al. (1995), but the scatter is very large. 
Fig. 4 shows that ~98% of the stars behind the dusty 
screen and a screen absorption of Av(scr) ~ 12 mag can 
explain the color excess observed along the gas bar. The 
regions next to the bar can be explained by a fraction of 
stars behind the screen of ^97% and an extiction reduced 
to Av(scr) ~ 7.5 mag. 

If we adopt an intermediate model with fhid = 97.5% for 
both the regions along and outside the gas bar, then we ex- 
pect the following differences in terms of surface brightness 



(mag/arcsec^) between the two zones: A/i^; w 0.27 ± 0.2, 
AfiH « 0.65 ± 0.35 and AfiK « 0.4 ± 0.2. Observation- 
ally, we obtain: A^r « 0.37 ± 0.22, A^h « 0.30 ± 0.10 
and AfiK ~ 0.15 ± 0.04. There is a fair agreement be- 
tween observations and expected values that, however, is 
worse at longer wavelength. This probably indicates that 
our model is too simplistic and, in particular, that a frac- 
tion of star is probably mixed with the dusty screen and 
that the dusty screen is not completely homogeneous as 
assumed. 

The visual extinction along the bar is about 12 mag- 
nitudes. If we assume the Galactic gas-to-dust ratio 
(Nh = 2.2 X 10^^ Ay cm~^), then the column density 
is Nh « 2 X lO^^cm"^. Note, however, that for active 
galaxies the standard conversion factor very likely un- 
derestimates the real N//, as discussed in Maiolino et al. 
(1999b). 



3.2. Evidence for streaming motions 

The kinematics of the gas in the presence of a barred po- 
tential has been studied by several authors (Athanassoula 
1992, Finer et al. 1995, Wada & Habe 1995, Friedh & Benz 
1993). Either gravitational torques due to the offset of the 
gas lanes with respect to the bar axis, or shocks along the 
gaseous lane, cause the gas on the leading side of the bar 
to lose its angular momentum and, as a consequence, to 
flow toward the nucleus. These motions translate into a 
strong velocity gradient across the leading side of the bar 
that should be observed in the Doppler maps of the emis- 
sion lines. Such a velocity pattern has been observed in 
several large scale stellar bars (eg. Laine et al. 1999, Re- 
gan et al. 1997). We have used our integral field near IR 
data to trace the velocity field of the molecular gas with 
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Fig. 5. — Color image: velocity (Doppler) field of the molecular gas as traced by the H2 S(l)(l-0) emission line at 2.12 fim, obtained by 
means of the 3D data. Green line: one of the H— K contours indicating the location of the dusty bar. h) Model of the velocity field expected in 
the case of absence of a barred potential (i.e. axisymmetric rotation), c) Difference between the observed velocity field and the axisymmetric 
model in b. 
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the H2 S(l)(l-0) emission line at 2.12 fira. Fig. 5 gives 
the radial velocity of the molecular gas as traced by the 
centroid of the H2 line and relative to the systemic velocity 
of 438 km/s, derived from the 21 cm line^. The color bar 
on the right-hand side of Fig. 5 translates the color coding 
into the radial velocity corresponding to the shift of the 
line. The green curve is one of the contours of the H-K 
map indicating the location of the gas bar. We excluded 
from our velocity mapping most of the regions that were 
poorly sampled by the dithering (i.e. with low effective 
integration time) since the signal to noise on the H2 line 
is too low in most of these regions. 

The most prominent features of this map are the strong 
velocity gradient and the highly redshifted gas along the 
leading side of the bar, indicating gas streaming toward 
the nucleus in agreement with the model predictions. This 
velocity field not only supports the reality of the gaseous 
bar, but directly shows inflow motions that probably con- 
tribute to the feeding of the active nucleus. 

The gas velocity pattern is similar the what observed in 
large scale stellar bars. Here the main difference, besides 
the much reduced size, is that the mass that exerts the 
torques generating the gas inflow is not dominated by the 
stars but by the gas. 

To show that the redshifted gas along the leading side is 
not due to noise or statistical fluctiations. Fig. 6 shows the 
spectrum of the H2 line along the leading side and com- 
pares it with the spectrum of the nucleus and with the 
spectrum of the central part of the bar (see figure caption 
for details). The line peak was normalized to one. The 
signal-to-noise on the line is high, and the shift of the line 
along the leading side is cleary systemic and not due to 
noise fluctuations. 

The streaming motions become obvious when compared 
to the rotational velocity field expected to dominate the 
kinematics in the absence of a barred potential. We have 



modeled the velocity pattern expected if the gravitational 
field is dominated by the stars that are traced by the K- 
band light (corrected for extinction), characterized by an 
average mass-to-light ratio M/Lk « 2.5 Mq/Lq (Maiolino 
et al. 1998). The resulting velocity field is shown in 
Fig. 5b. In the outer regions (> 100 pc, where the bar 
should not play a role) , it is in agreement with the observed 
velocity field (Maiolino et al. 1998, Storchi-Bergman et al. 
1999). Fig. 5c shows the deviations of the observed ve- 
locity field relative to the rotational field expected in the 
axisymmetric case: the most obvious deviation is the strip 
of redshifted gas along the leading side of the gaseous bar, 
in agreement with model predictions and as observed in 
large scale stellar bars. 

The velocity maps (both Figs. 5 a and c) also show a 
strong gradient along the dusty spiral arm that extends 
from the southern end of the bar to the east, also as pre- 
dicted by models (see eg. Finer et al. 1995). There is 
also a region of redshifted gas to the east of the nucleus; 
some numerical simulations show the formation of a ve- 
locity gradient in this region (Athanassoula 1992, Finer et 
al. 1995): indeed this region is located about at one of 
the aphelia of the bar, where a significant population of 
X2 orbits can reduce the average velocity of the gas with 
respect to the circular axisymmetric motions. 

The largest velocities of the streaming gas along the bar 
are expected to occur close to the center (near the peri- 
helia), as is actually observed in large scale stellar bars. 
Such high velocity gas is also observed in the nuclear re- 
gion of Circinus. In Fig. 7a the green contour is a channel 
map of the H2 emission with velocity Doppler shift be- 
tween -1-200 and -1-260 km/s, while the color image is the 
[SiVI] map that will be discussed in Sect. 0. There is clear 
evidence for highly redshifted gas at the end of the gaseous 
bar, close to the nucleus. In this region the deprojected 
velocity is about 300 km/s. 



^The nuclear velocity of H2, as well as of other low ionization lines and of the stellar features (Fig. 6 of Maiolino et al. 1998), shows a small 
redshift of about 20 km/s. This redshift was also observed in high resolution echelle spectroscopic observations (Oliva, private communication). 
This redshift might arise from the definition of the velocity for the 21 cm line. Otherwise, it might reflect small, systemic differences between 
the central region and the outer parts of the galaxy. 
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The velocity field of the central region (20-40 pc) is 
however dominated by a normal rotation pattern: here 
the isovelocity lines are parallel to the minor axis. This 
nuclear rotation might be associated with the disk-like fea- 
ture observed in the H-K color map. 

Summarizing, the molecular gas velocity field shows 
clear non-circular motions. The general behavior of the 
kinematics can be ascribed to streaming motions due to 
the gaseous bar. Obviously, we cannot rule out more com- 
plex situations, especially concerning irregularities and 
details of the velocity field. Higher angular and spectral 
resolution and a more detailed modelling of the gas kine- 
matics in this region are required to further support our 
findings and/or provide new insights on the details of the 
gas dynamics. 



4. THE INNERMOST REGION OF THE IONIZATION CONES 
4.1. The 10 pc scale 

The gas inflow through the bar not only brings mate- 
rial close to the nuclear black hole, but probably plays 
an important role also in other phenomena of the nuclear 
activity such as the ionization cones and the star forma- 
tion. Fig. 7a shows the NICMOS map of the [SiVI] coronal 
emission line. The [SiVI] emission extends in the ioniza- 
tion cone traced by the [OIII] optical images (Marconi et 
al. 1994a), as indicated by the dashed lines. The presence 
of [SiVI] emission also on the eastern side of the nucleus 
provides evidence for a counter-cone, whose existence was 
also inferred from the larger scale radio maps (Elmouttic ct 
al. 1995, 1998c). The foreground extinction due to the nu- 
clear dusty feature (disk?) is very likely responsible for the 
reduced [SiVI] brightness of the counter cone. If the nu- 
clear dusty feature produces an extinction of Ay ^ 9 mag 
(as inferred from Fig. lb), this would imply Argjvi] ^1.1, 
fully accounting for the brightness difference between the 
two cones, about a factor of 3. 

The [SiVI] emission is not azimuthally uniform within 
the (NW) ionization cone; the southern edge appears much 
brighter than the rest. It is likely that the gas trans- 
ported by the bar enters the ionization cone in this region. 
This interpretation is supported by the distribution of the 
molecular gas shown in Figs. 7 a and 6, where we report 
the 3D maps of the integrated Il2S(l)(l-0) emission line 
and the highest velocity component of the gas traced by 
this line. These maps show the presence of large quantities 
of gas just south of the ionization cone, that is flowing at 
high velocities towards the southern edge of the ionization 
cone. Next to the edge of the cone the H2 emission is high- 
est because here the molecular gas is more exposed to the 
nuclear radiation, but also because the gas density is high 
and the II2 has not had time to dissociate fully. Indeed, 
the R-H color map (Fig. 3) clearly shows a prominent dust 
lane extending along the southern edge of the ionization 
cone, coincident with the H2 emission. As soon as these 
dense clouds enter the ionization cone they are ionized and 
emit the [SiVI] line, but also Paa (Fig. 7c); this gas is very 
dense and therefore the ionized gas emission is very high. 
Another dust lane in Fig. 3 traces the northern edge of the 
ionization cone. 

Another possible explanation of the emission line and 
dust lane morphology is that these features are simply the 



effect of the orientation of the ionization cone with re- 
spect the nuclear part of the galactic gaseous disk, which 
might be warped on these small scales (Quillen et al. 1999, 
Mulchaey et al. 1996, Schreier et al. 1998). The brighter 
part of the ionization cone might correspond to the region 
closer to denser gas of such nuclear disk. The only diffi- 
culty with this model is to obtain a projected velocity in 
excess of 200 km/s, as observed in Fig. 7a. 

As shown in Figs. 7 a and b the molecular clouds en- 
ter the ionization cone very close {'^ 10 pc) to the nu- 
cleus, where the radiation pressure is very high. As dis- 
cussed in Binette (1998) the radiative acceleration is given 
by a^ad « (n/1200cm-3)(Uo/0.125)(3.1 x 10"^ + 5.8 x 
10~^^d) [cm s~^], where n is the gas numerical density, 
Uo is the ionization parameter and /iD the dust content 
with respect to the solar value. The density of the coronal 
line region was estimated to be about 5000 cm~^ based 
on the infrared [NeV] doublet (Moorwood et al. 1996). 
Oliva et al. (1999) determined Uq for a single cloud in the 
ionization cone; we scaled this value to the distance of 10 
pc and to the density of 5000 cm~^, i.e. Uq « 1. With 
these values arad ~ 3 x 10~^cm s~^ for a dust content of 
fiD = 0.1. Maiolino et al. (1998) set an upper limit to the 
mass of a nuclear black hole in Circinus of 4 x lO^M©. The 
gravitational pull at a distance of 10 pc is estimated to be 
4 X 10~^cm s~^, i.e. about one order of magnitude lower 
than the radiative acceleration. Therefore, these clouds 
just driven into the ionization cone will be pushed out- 
ward to form an outflowing component of the Coronal Line 
Region and of the NLR. This outward acceleration of the 
coronal line clouds, as well as the tendency of the clouds 
to be ablated under the strong radiation pressure (Math- 
ews 1986), is probably responsible for the lower emissivity 
of the [SiVI] line in the northern part of the cone. The 
denser clouds injected into the cone will survive ablation 
and will be accelerated like "bullets". Veilleux & Bland- 
Hawthorn (1997) discovered the presence of supersonically 
ejected clouds (Vojoct > 100 km s^^) in the Circinus NLR, 
which they ascribe to some explosive nuclear event. We 
claim that such ejected "bullets" are the result of the radia- 
tive acceleration of the dense clouds that have been driven 
into the ionization cone by the gaseous bar or by the nu- 
clear disk. As suggested by Binette (1998), the radiation 
pressure-driven outflow would explain the blueshift of the 
coronal lines observed in Circinus (Oliva et al. 1994) and 
in other Seyfert galaxies (Pension et al. 1984), as well as 
the dependence of the blueward asymmetry of the narrow 
lines on the ionization potential and on critical density 
(DeRobertis & Shaw 1990, Dahari & De Robertis 1988, 
Whittle 1992, Heckman et al. 1981). A potential problem 
of the "outflow model" would be that new gas must be con- 
tinuously supplied to the inner part of the cone to replace 
the outflowing gas. The gas flow provided by the gaseous 
bar or by the nuclear disk could to be the mechanism re- 
sponsible for feeding the inner part of the ionization cone 
(a more quantitative discussion is given in sect. 6.4). 

On the other hand, a large fraction of the gas in the 
outer regions of the NLR (100 pc scale), characterized 
by lower ionization emission species, is probably part of 
the normal galaxy disk illuminated by the nuclear source. 
Indeed, the kinematics of the outer parts of the ionization 
cone are also characterized by a rotational component (El- 
mouttie et al. 1998b, Veilleux & Bland-Hawthorn 1997). 
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Fig. 7. — a) NICMOS map of the [SiVI]1.96 /^m line (color) along with the 3D channel map of the H2 2.12 /^m line with the highest 
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Paa emission line, d) NICMOS map of the [Fell] 1.64/xm emission line. The green dashed lines indicate the orientation of the light cone 
observed in [OIII]. 



4.2. The parsec scale 

Another interesting feature of the emission hne maps in 
Fig. 7 is that both the [SiVI] and [Fell] hues peak on the 
nucleus: although both show extended emission, a large 
fraction of the integrated flux (^ 70% for [SiVI]) comes 
from an unresolved source coincident with the 2.2/im con- 
tinuum peak, i.e. the location of the active nucleus. The 
shift of the [SiVI] peak observed in the lower resolution 
data of Maiohno et al. (1998) was due to a blend of the 
nuclear component with the extended component to the 
west. Although, a fraction of the coronal lines might be 
excited by shocks, most of the coronal line emission must 
be ascribed to photoionization (Oliva 1997). However, ac- 
cording to photoionization models, the [SiVI] is expected 
to peak at the outer edge of the He"*"^ Stromgren sphere, 
while the [Fell] should peak right after the H+ Stromgren 
sphere. With the luminosity of Circinus and a density of 
10^ cm~'^ (i.e. even higher than what is estimated through 
the [NeV] infrared lines) Oliva et al. (1999) derived that 
the [SiVI] should peak at about 30 pc from the nucleus. 
The [Fell] emission should peak even further out. There 
are a few possible explanations of this inconsistency: 



1 



spheres. This might be the case for the [Fell] emis- 
sion. Yet the clouds responsible for the emission of 
most coronal lines (including [SiVI]) appear charac- 
terized by little dust content and, therefore, dust 
probably plays a minor role in the ionization struc- 
ture of these clouds. 

2. Another possibility is that the gas that emits most 
of the coronal lines is significantly denser, shrink- 
ing the size of the He^^ sphere. More specifi- 



cally R(He+2) oc N, 



-2/3 



Reducing the radius of 



As suggested by Oliva et al. (1999), the ionized 
clouds are dusty. The ionization parameter in these 
regions is very high (U'^ 1) and, therefore, the dust 
competes with the gas in absorbing UV photons 
thus reducing significantly the size of the Stromgren 

'^Oliva et al. show that very likely most of the Fe is in the gaseous form in the coronal line region. 



the [SiVI] emitting shell from 30 pc to less than 
10 pc (which would make the emission unresolved) 
would require an electron density larger by a fac- 
tor of 5, i.e. No = 5 X lO'^cm^'^. However, this 
density is one order of magnitude higher than what 
measured from the [NeV] doublet (Moorwood et al. 
1996). Also, the thickness of the [SiVI] emitting 
shell (Marconi et al. 1994b), relative to R(He+2), 
is inversely proportional to Ng and, therefore, the 
[SiVI] line would be significantly fainter than other 
coronal lines that form inside the He"*"^ sphere, i.e. 
higher ionization species. This line of reasoning ap- 
plies as long as the ionization structure has a quasi- 
spherical geometry (Oliva 1997). If a nuclear popu- 
lation of very dense clouds exists, more specifically 
with Ne > lO^cm"^ (but still lower than the critical 
density of ~ 10®cm~^), then the geometry of the ion- 
ization structure is quasi-plane parallel (R > 10 AR) 
and the region of [SiVI] emission extends over the 
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whole Hc+^ Stromgrcn region. This would solve the 
problem of the ratio between [SiVI] and the higher 
ionization coronal lines, but the lower density mea- 
sured with the [NeV] lines remains an issue. 

3. Increasing the filling factor would also reduce the size 
of the emission region (R(He'*'^) oc f^^/^), but this 
would also affect the coronal line ratios similarly to 
the change in density discussed above. 

4. A final, interesting possibility is that the nuclear UV 
emission is intrinsically non-isotropic (Oliva et al. 
1999), as is expected to be the case if the emission 
arises from an optically thick accretion disk. In the 
latter case the azimuthal clouds are illuminated by a 
flux that is significantly lower than the face-on flux 
(Lagn oc {cos 9)-^, Laor & Netzer 1989). The He+2 
Stromgren radius of the azimuthally distributed gas 
would be much smaller than that estimated above 
and could very well appear point like at our resolu- 
tion. 

5. NUCLEAR STAR FORMATION 

The Paa and the [Fell] emission appear much more dif- 
fuse and extended® than the [SiVI]. Both emission lines 
extend well outside the ionization cone and preferentially 
along the major axis of the galaxy, although the Paa also 
has a component in the ionization cone similar to that ob- 
served in the [SiVI] map. Very likely the diffuse/extended 
emission traces moderate star forming activity occurring 
within the central 50 pc as a result of the large amount of 
molecular gas driven into this region by the gaseous bar. 
The Paa emission traces the UV radiation from young 
OB stars. Maiolino et al. (1998) used the circumnuclear 
Br7 emission and the mass-to-light ratio to model the star 
formation history, and infer that the nuclear stellar popu- 
lation has an age of about 10* yr. 

Fig. 7c shows that the star formation history of the nu- 
clear region is probably more complex. The Paa map 
shows the presence of knots and non-axisymmetric emis- 
sion: a chain of three HII region is present at about 20 
pc to the north-east of the nucleus, while a burst of star 
formation is present to the south-west, right at the nuclear 
end of the leading side of the gaseous bar. These are the 
locations of the aphelia of the bar, where orbit crowding is 
indeed expected to enhance the star formation (Wada and 
Habe 1992). There is no evidence for star forming activity 
along the bar. Although large concentrations of gas are 
expected along the leading side of the bar, star formation 
might be inhibited by the strong shear forces (Elmegreen 
1979, Kenney & Lord 1991). 

The diffuse [Fell] emission probably traces gas shocked 
by supernova explosions. A bright [Fell] knot, about 30 pc 
to the north of the nucleus very likely traces a single young 
supernova remnant. The [Fell] knot is unresolved (i.e. R< 
4 pc), thus implying that the SNR must be younger than 
about 100 yr. There is also a prominent [Fell] extension 
to the south of the nucleus, slightly offset with respect to 
the Paa emission; shocks due to the gas streaming from 
the gaseous bar might also play a role in this region. 

6. DISCUSSION 



6.1. Plausibility of the gas bar and open issues 

Although the presence of a nuclear gas bar appears con- 
vincing both from the morphological and the kinematical 
points of view, we can check whether the physical con- 
ditions of the gas meet the requirements necessary to de- 
velop the gravitational instabilities that should form a bar. 
According to models, a gaseous bar should form when the 
mass of the gaseous nuclear disk is larger than 10 (Wada & 
Habe 1992) to 20% (Shlosman et al. 1989), of the dynam- 
ical mass. C0(2-l) millimetric data (Aalto et al. 1995, 
Johansson et al. 1991) indicate that the gas mass en- 
closed within the central 400 pc is about 6 x IO^Mq, if the 
Galactic conversion factor is assumed (Kenney & Young 
1989). Based on the kinematical information provided by 
Maiolino et al. (1998) we estimate that the dynamical 
mass inside the same region is about 1.2 x 1O^M0, i.e. 
only a factor of two higher than the mass in molecular 
gas. Therefore, the physical conditions of the gas in the 
nuclear region are indeed appropriate for the formation of 
a gaseous bar (although this argument is weakened if the 
conversion from CO to gas mass is less than the standard 
Galactic value, as sometimes appears to be the case in 
extreme starbursts). 

In the H-K color map (Fig. lb) only one side of the 
gaseous bar is observed; the north-west counterpart is not 
seen. One possible explanation is that the other side of 
the gaseous bar is present but it is not observed in the ex- 
tinction map because the contrast of absorption features 
is much lower on the far side of the disk (more foreground 
stars and less stars behind). Otherwise the north-west side 
of the bar might be intrinsically weaker or even absent. 
Since gaseous bars are formed via gravitational instabil- 
ities, they are expected to be highly inhomogeneous and 
possibly asymmetric, especially on such small scales. 

Another important issue to address is whether these 
gaseous bars are rare cases or they could be a fuelling 
mechanism common to most active galaxies as claimed by 
Shlosman et al. (1989). If the latter is the case then 
the question is why such nuclear gaseous bars were de- 
tected only in a few active nuclei (Mulchaey & Regan 1999, 
Hcnkel et al. 1991, Ishizuki et al. 1990), or none if we 
only consider nuclear bars as small as the one in Circinus. 
These gaseous bars can be identified from molecular gas 
emission (eg. CO millimetric maps) and/or through dust 
absorption traced by the color maps, as we have done in 
this work. Circinus is closest Seyfert 2 known and the pro- 
jected angular size of the gaseous bar is 2 arcsec. Most of 
the other Seyfert galaxies are further away and, therefore, 
a nuclear gaseous bar with similar size would have a much 
smaller angular size, which would make it hard to detect. 
The average distance of the Seyfert galaxies in the sample 
of Maiolino & Rieke (1995) is 32 Mpc and the projected 
size of the nuclear gaseous bar would be a few tenths of arc- 
second, unresolved even at the NICMOS resolution (at 2 
fim) and at the resolution of the current millimetric inter- 
ferometers. Also, we emphasize that the color map tech- 
nique to search for gaseous bars works as long as infrared 
filters are used. In these heavily obscured nuclei the opti- 
cal images sample only the outer (foreground) parts of the 
nuclear region and therefore do not provide the required 
information even if compared with an infrared image to 



^The [Fell] morphology observed in our NICMOS image is consistent with the lower resolution map obtained by Davies et al. (1999). 
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obtain a color map. Indeed, as shown in Fig. 3 and as 
discussed in sect.3.1, the R-H color map of Circinus does 
not show evidence for the dusty bar detected in the H-K 
color map. This result may partly explain the shortage of 
gaseous bars in Regan & Muchaey (1999), a work which 
included images taken through the above filters of a sam- 
ple of Seyferts. Summarizing, whether nuclear gas bars 
are a common feature in Seyfert galaxies or not remains 
an open question, since such structures might have been 
missed in previous studies both because of limited angular 
resolution and/or inappropriate techniques. 

6.2. Fuelling the active nucleus 

As discussed in Athanassoula (1992) and Regan et al. 
(1997), probably only a small fraction of the streaming 
gas does actually result in an inflow (10%-30%). We can 
estimate the streaming flux of the gas along the leading 
edge of the nuclear bar in Circinus and compare it with 
the inflow requirements to feed the active nucleus. The 
gas streaming flow along the leading edge of the bar (an 
upper limit to the real inflow rate) is given by: 



M« 



Nh W Vstroam mp 



(3) 



where w is the width of the streaming gas lane {^ 20 pc) , 
mp is the proton mass, Vstream is the streaming gas veloc- 
ity along the leading edge of the gaseous bar and Nh the 
gas colu mn density in the same region ('--^ 2 x lO^^cm"^, 
sect.3.1). To determine Vstroam we have to correct the ob- 
served redshifted velocity for various projection effects and 
for the contribution due to the pattern circular motion, as 
illustrated in Fig. 2. More specifically: 



Vs. 



VD(obs)/ cos i - Vpattern siu I 
cos Of 



(4) 



where Vd is the Doppler velocity shift observed along the 
leading edge of the bar, is the deprojected angle between 
the galaxy minor axis^ and the line connecting the nucleus 
to the specific point along the leading edge of the bar where 
the redshifted gas is observed (see Fig. 2), Vpattorn is the 
velocity component due to pattern speed at t he same lo- 
cation, while i and a were defined in Sect. |3.1| . The outer 
edge of the gaseous bar is located at about 160 pc from 
the nucleus, probably the location of the corotation radius. 
At this distance the rotational velocity is 150 km/s (depro- 
jected for inclination, MaioHno et al. 1998), thus giving 
a pattern speed of ripattern = 0.92 km s~^ pc~^. If we 
consider the streaming gas located at about the mid-point 
of the gaseous bar, then the deprojected distance from the 
nucleus is about 90 pc, hence Vpattorn = 83 km/s, and 
(j) = 25°. The contribution from the circular pattern ve- 
locity (Vpattorn siu 0) is very low, as expected since the bar 
is oriented almost along the line of sight. With these val- 
ues we estimate Vstroam = 83 km/s, hence from Eq. p3 we 
derive a streaming flow rate of Mstroam ~ 0.25 M0yr~^. 

The luminosity of the active nucleus is about 2 x 10^° Lq, 
or lower if the nuclear emission is intrinsically anisotropic 
(Oliva et al. 1999), thus implying an accretion rate of 

Magn = LAGN/c^e w lO-^Meyr-i (5) 

®More specifically the direction orthogonal to the line of nodes and 



(where we assumed a radiation efficiency e = 0.1). The 
streaming flow estimated above is 25 times higher than 
the AGN accretion rate. Therefore even if only 10% of 
the streaming gas actually flows into the nucleus, this will 
provide 2.5 times as much fuel as required to power the 
activity. 

6.3. Fuelling the nuclear star formation 

As discussed above, Paa emission traces regions of re- 
cent star formation, while the smoother near-IR contin- 
uum emission traces the older stellar population evolved 
from past bursts. We can estimate the current star for- 
mation rate from the Paa emission due to stellar activity 
integrated in the nuclear region. We subtracted the con- 
tribution to the Paa emission due to the active nucleus 
by normalizing the [SiVI] map to the Paa on the nucleus; 
this is actually a conservative subtraction of the AGN Paa 
flux, since a fraction of the unresolved nuclear Paa emis- 
sion, as well as part of the Paa emission within the cone, 
might be due to star formation as well. The flux of the 
AGN-subtracted Paa integrated within a radius of 30 pc 
(= 1".5) is 1.2 xlO^^'^ ergs~^ cm~^. If we assume that the 
extinction derived from the H-K color map in Fig. 1 also 
applies to the nebular emission, then on average Ay ~ 9 
mag, hence Apaa ~ 1.3 mag. Therefore the extinction 
corrected Paa flux is 4 x 10"^"^ erg s~^ cm~^. At the dis- 
tance of 4 Mpc, this flux translates into a ionizing photon 
flux Qi = 1.4 X 10^^ s^^. We employed the stellar popula- 
tion synthesis code described in Sternberg & Kovo (1999) 
to derive the star formation rate required to produce this 
ionizing photon flux. By adopting a Salpeter initial mass 
function and an upper mass cut-off of 60 M© we obtain a 
star formation rate of ^^ 10^^ Moyr^^. This star forma- 
tion rate is relatively low and similar to the accretion rate 
of the active nucleus. T he a mount of gas inflow provided 
by the gaseous bar (sect.|]^) can feed both the active nu- 
cleus and the star forming activity in the nuclear region. 

6.4. Fuelling the ionization cone 



As discussed in sect. 4.1 



part of the gas observed in the 
ionization cone is in circular motion in the galactic plane, 
while a fraction of the gas is outflowing, probably as a 
consequence of the nuclear radiation pressure. We spec- 
ulated that the inflow due to the gas bar might supply 
to the nucleus the material that is being ejected through 
the ionization cone. Here we address this possibility more 
quantitatively. Veilleux & Bland-Hawtorn (1998) estimate 
the kinetic energy in the outflowing gas and in the ejected 
filaments to be about K «:! 10^^(Nc/100cm~^) erg. By ap- 
proximating with a spherical outflow, the kinetic energy 
can be converted into a mass outflow rate according to the 
formula Moutfl ~ 2K/(VoutflR), where Voutfl is the (depro- 
jected) outflow velocity and R is the average distance of 
the outflowing gas detected in the images. By assuming 
Voutfl « 200 km/s (Veflleux & Bland-Hawtorn 1998) at 
an average distance of about 100 pc, and by assuming an 
average density No « 500cm^^ (Oliva et al. 1999), we 
obtain Moutfl ~ lO~^M0yr~-'^. Although this estimated is 
very rough, it shows that the outflowing component of the 
ionization cone requires some degree of fuelling, and that 
the inflow due to the gas bar can account for this fuelling. 

on the plane of the galaxy. 
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7. CONCLUSIONS 

We present NICMOS/HST images and integral field 
near-IR spectroscopy of the nuclear region of the closest 
Seyfert 2 galaxy Circinus. The H-K color map shows a nu- 
clear dust lane that extends radially. We propose that this 
dust lane traces a nuclear gas bar that might be responsi- 
ble for feeding the active nucleus. The velocity field traced 
by the H2 line at 2.12 /im shows enhanced radial motions 
along the leading edge of the gas bar, in agreement with 
model predictions. The latter finding not only supports 
the identification of the gaseous bar, but directly shows 
gas inflow motions that probably fuel the AGN. The es- 
timated streaming flow along the leading edge of the bar 
is about 0.25 Moyr~^. Probably no more than 10-30% 
of this flow does actually result in an inflow rate, that is 
however enough to fuel the active nucleus (whose accre- 
tion rate is about 10~^ MQyr~^). We cannot rule out 
more complex situations to explain details and irregular- 
ities of the velocity field; higher resolution spectroscopic 
data are required to further support our findings and/or 
provide new insights on the details of the gas dynamics. 

Nuclear gaseous bars similar to the one in Circinus might 
be a fuelling mechanism common to many active nuclei 
hosted in gas rich galaxies. Possibly, limited angular reso- 
lution and/or inappropriate techniques might be responsi- 
ble for the shortage of nuclear gas bar detections in other 
active nuclei. 

We also used the NICMOS data to study the innermost 
part of the Narrow Line Region and the Coronal Line Re- 
gion. The coronal line [Si VI] 1.97 /im reveals the presence 
of a counter-cone that was not observed at optical wave- 
lengths, very likely because of obscuration. The [SiVI] and 
Paa emission is not distributed uniformly within the ion- 
ization cone but extends preferentially along the southern 
edge of the cone. This enhanced emission seems associ- 
ated with dense gas just outside the cone, traced by the 



H2 (molecular gas emission) and R-H (dust absorption) 
maps. This gas is flowing towards the cone at velocities in 
excess of 200 km/s. We speculate that the gas streaming 
from the gas bar enters the ionization cone, contributing 
significantly to the Coronal Line Region and to the inner- 
most part of the Narrow Line Region. This dense gas is 
probably accelerated by the radiation pressure producing 
the outfiowing gas and the ejected clouds observed in the 
Narrow Line Region. 

Besides the extended component, both the [SiVI] 1.97 
/im and the [Fell] 1.64 /im emission are characterized by a 
nuclear intense unresolved component, while photoioniza- 
tion models would expect these lines to peak at a few 10 pc 
from the nucleus. This finding poses important constraints 
on the photoionization models. 

Both Paa and [Fell] emission extend over the central 60 
pc and also outside the ionization cone. Most probably 
this extended emission traces recent star formation in the 
nuclear region. The Paa extended emission is irregular 
and "knotty" , different than the smoother near-IR contin- 
uum. This Paa emission probably traces recently formed 
HII regions, therefore pointing to a complex nuclear star 
formation history. We estimate an average star formation 
rate in the nuclear 60 pc of ~ 10~^ MQyr"-'^, which could 
be sustained by the inflow rate doe to the gaseous bar. 
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